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Abstract 

We suggest searching for excited charmonium states in e + e~ annihilation via double charmonium 
production at ^/s = 10.6 GeV with B factories, based on a more complete leading order calculation 
including both QCD and QED contributions for various processes. In particular, for the C=+ 
states, the \co{nP) (n=2,3) and r] c (mS) (m=3,4) may have appreciable potentials to be observed; 
while for the C=- states, the r] c h c production and especially the Xcih c production might provide 
opportunities for observing the h c with higher statistics in the future. A brief discussion for the 
X(3940) observed in the double charmonium production is included. 

PACS numbers: 12.40.Nn, 13.85.Ni, 14.40.Gx 
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I. INTRODUCTION 



Charmonium spectroscopy has become a challenging topic in hadron physics and QCD, 
because of the recent findings of possible new charmonium states (for recent experimen- 
tal and theoretical reviews and discussions, see e.g. and references therein). 
Among others, for the puzzling state X(3872), possible assignments of e.g. the 2 3 Pi and 
l 2 Di charmonium states and the charm-molecule have been suggested (see, e.g. j^J for a 
comprehensive review), and it will be helpful to search for those states in other experiments 
and to clarify these assignments; The measured mass splitting between ip(2S) and r] c (2S) 
is about 50 MeV, which is smaller than some theoretical predictions, and it is certainly 
useful to search for the r] c (3S) to see what will be the mass splitting between the ip(3S), 
which could be the observed t/>(4040), and the r] c (3S). This may be particularly interesting 
since according to some potential model calculations the i] c (3S) could lie above 4040 MeV 
(see, e.g. in [6| the mass of r] c (3S) is predicted to be 4060 MeV). And the tp(3S) mass 
could actually be lowered by coupling to the nearby D*D* decay channels (note that the 
energy level spacing between ip(3S) = ?/>(4040) and ip{2S) = ?/>(3686) is smaller than that 
between ip(4:S) = ip(4414) and i/j(3S) = ^(4040), which is in contradiction with potential 
model calculations unless the coupled channel effects are considered or the assignments for 
ip(3S) = ?/>(4040) and ip(4S) = ^(4414) are incorrect. The mass spectrum of excited char- 
monium states will certainly provide important information on interquark forces and color 
confinement. In addition, studies of the decay and production of these states will also be 
very important in understanding the underlying theory of strong interaction -perturbative 
and nonperturbative QCD in view of many outstanding puzzles in charmonium physics. 

B meson decays have proven to be very useful processes to find new charmonium states. 
Aside from the B meson decay, e + e~ annihilation at = 10.6 GeV could also be a 
very usefulprocess in finding the excited charmonium states, since the recent Belle exper- 
iments [3, Is} have found unusually strong signals for the double charmonuim production 
from the e + e" continuum, e.g., e + e _ — > J/4>r] c , J/ipXco, J/4>Vc(2S) and e + e~ — > ip(2S)rj c , 
ip(2S)xco, i/j(2S)r) c (2S) . Theoretically, the calculated cross sections for these processes based 



on the leading order Non-Relativistic QCD(NRQCD) (or more generally pertur 
(pQCD)) are about an order of magnitude smaller than the experiments 9l. llOl. 



Dative QCD 



This is 



a big issue in charmonium physics and NRQCD, and it still remains to be further clarified 
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though many considerations are suggested to understand the large production rates in both 
exclusive and inclusive charmonium production via double charm pairs in e + e~ annihilation 



(the theoretical 



color-singlet 



14 



predictions for the inclusive J/^cc production cross section with the 



lq as well as color-octet [1 61] contributions are also much smaller than 



the Belle data). Despite of these disagreements, however, we find that the calculated rela- 
tive rates of the double charmonium production processes are roughly compatible with the 
Belle data (e.g. the production cross sections of i] c , r) c (2S), and Xco associated with J/ijj 
and i/j(2S) are much larger than that of Xd an d Xa)- So, we may use the same method as 
in our previous work to calculate the production rates for the excited charmonium states 
in e + e~ annihilation into double charmonia, but mainly pay attention to the relative rates 
for these production processes. We hope the calculation will make sense in predicting the 
relative production rates for those excited charmonium states, and can be tested by experi- 
ments. This will be useful not only in the search for those excited charmonium states, but 
also in understanding the production mechanism itself. If the predicted relative production 
rates turn out to be consistent with experiments, it is likely that the NRQCD factorization 
treatment for these processes probably still makes sense and only an overall enhancement 
factor is needed and should be clarified in further theoretical considerations (including QCD 
radiative corrections, relativistic corrections, and other nonperturbative QCD effects). In 
the last section we will have a discussion on recent developments in this regard. In the 
following, we will calculate the leading order production cross sections for various excited 

n 

charmonium states in e + e annihilation at y/s = 10.6 GeV in the same way as in [1 0f | . 
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II. FORMULAS AND CALCULATIONS 



Following the NRQCD factorization formalism 171]. the scattering amplitude of double 
charmonia production can be described as 



x 



A(a + b^QQC s ^L J ^(p 3 ) + QQC s ^L J )(p A )) = Jc^y/cl £ EEEEE 

L^ z S^z siS2 jk L Z S Z U 

X (s 

x (s 3 ;s 4 I SS Z )(LL Z ;SS Z \ JJ z )(3l;3i | 1) 

X(a + b - Q,-(f ) + Q fe (f ) + Q,(a) + Q.(B)) (L = 5), 

e* (L z )^(a + b - g,(f ) + Q fc (f ) + Q,(f ) + Qi(f )) = P), 

(1) 

where (3j;3k | 1) = 5 jk /y/W c , (3Z; 3i | 1) = 5u/\JW c , (si;s 2 I £V>SW , («3; s 4 I 
S'S'z) , (L^L^z] S^pS^ z I J-^J^z) and (LL Z ;SS Z \ JJ Z ) are respectively the color-SU(3), 
spin-SU(2), and angular momentum Clebsch-Gordon coefficients for Q<5 pairs projecting 
out appropriate bound states. A(a + b — ► + Ofc(^f) + Qi(^) + Oi(^f)) is the scat- 

tering amplitude for double production and A a is the derivative of the amplitude with 
respect to the relative momentum between the quark and anti-quark in the bound state. 
The coefficients Cx . and Cl can be related to the radial wave function of the bound states 
or its derivative with respect to the relative spacing as 

^ = ^(0)| 2 , C p = ^\R' P (0)\ 2 - (2) 
We introduce the spin projection operators Pss z {P-,q) as 

HQ 

p ss z (p,q) = ^2(si;s 2 \SS z }v(^ - + q; s 2 ). (3) 

Expanding P$s z {P, q) in terms of the relative momentum q, we get the projection operators 
and their derivatives, which will be used in our calculation, as follows 

PisAv,®) = ^4*{S z ){£> + 2m c ), (4) 
Poo(p,0) = ^= l5 ( J p+2m c ), (5) 
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rom a potential model calculation (see the QCD (BT(Buchmuller- 
: |i? ls (0) | 2 = 0.810GeV 3 , \R 2S (0) | 2 = 0.529GeV 3 , \R 3S { 



P&(p,0) = —^-^jA (6) 
2y 2m c 

P 1 %(P^) = ^^h a f(S,)U + 2m e )-(f-2m e ) t(S z )j a }. (7) 

We then get the following expressions and numerical results for various processes of double 
charmonium production in e + e~ annihilation at y's = 10.6 GeV. In the calculation of the 
short distance coefficients, the quark and anti-quark are all on mass shell, and the meson 
masses are taken to be m 3 = m 4 = 2m c . The input parameters are \/s = 10.6GeV, m c = 
1.5GeV, a s = 0.26 (corresponding to /i = 2m c and A^,=338 MeV), and the wave functions 
at the origin are taken 
Tye)) model in Ref. 

0.455GeV 3 , |i?' 1P (0)| 2 = 0.075GeV 5 , |^ P (0)| 2 = 0.102GeV 5 , and 1^(0) | 2 = 0.015GeV 7 



A. The tp(nS)rj c (mS) production 

In the two S-wave (nS and mS) case, the cross section for e + + — * 7* — > H\{nS) + 
H 2 {mS) is given by 

( x(e + (p 1 ) + e-(p 2 ) -> H x {p z ) + H 2 {p A )) = 

2iraW s \R ns (0)\ 2 \R m s(0)\Ws 2 ~ 2s(mj + mj) + (m 2 f 1 ,,-v |2 , fl (R , 

L m dcos9 > (8) 

where 6 is the scattering angle between p~[ and p~s, \M\ 2 is as follows 

2 16384m c 2 (t 2 + n 2 - 32mj) 

|M ^' " (8m*-t-u)s • (9) 

Here the Mandelstam variables are defined as 



m 



S=(P1+P2)\ (10) 

t = (ps - Pi y = + + ^ 2 -2s(m 2 + m 2 ) + K-m 2 ) 2 ^ ^ (n) 

The cross sections for the double S-wave charmonium production are listed as follows 

a(e + + e- -> J/ip(lS) + Vc (lS) [r/ c (25), ^(35)]) = 5.5[3.6,3.1] fb, (13) 



a(e + + e~ -> ^(2S) + ^(15) [^(23), 7? c (3£)]) = 3.6[2.3, 2.0] fb, (14) 

where the values in the brackets are the cross sections respectively for r] c (2S) and r] c (3S) 
with recoiling -0 mesons. 

B. The ij}{nS)Xcj{mP){J = 0, 1, 2) production 

In the one spin-triplet S-wave (nS) and one spin-triplet P-wave (mP) case, the cross 
section for e + + e~ — > 7* — ► H^nS) + H 2 {mP) process reads 



<r(e + (pi) + e-(p 2 ) -> ^(pg) + #2^4)) 

-I- m 2 ) 4- (m 2 — m 2 ") 2 f 1 _ , 

|M| 2 dcosfl. (15) 



27ra 2 a s 2 |i?„ s (0)| 2 |^ p (0)| V« 2 - 2s(m§ + m 2 ) + (m 2 - m 2 ) 2 yi 



27m 2 s 2 
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Here for the production of spin-triplet states ipXcj, \M\ 2 is given in Eq. (Tl6|) for x c o, in 
Eq. flUDfor Xd and Eq. (TUD for Xc2 , 



|M too | 2 = 2048(90112m^ - 74752m^- 74752m^ + 23360m^ 2 + 43136m^ 



+ 23360m V - 3152m 4 t 3 - 7600m 4 t 2 w - 7600m 4 tw 2 - 3152m 4 w 3 



+ 162m 2 t 4 + 4Um 2 c t 3 u + 564m 2 tV + 444m 2 tw 3 + 162m 2 w 4 

3„.2 q+2„,3 + „.4\ //o„7 2\ 



- t 4 u-3t 3 u 2 -3t 2 u 3 -tu 4 )/{3s 7 m 2 c ). (16) 

|M tel | 2 = 32768(1792m* + 256m 6 c t + 256m[!u - 56m 4 c t 2 - 64m 4 tw - 56m 4 w 2 - 4m 2 t 3 

-20m 2 t 2 w - 20m 2 c tu 2 - 4m 2 w 3 + 1 4 + 2t 3 u + 2t 2 u 2 + 2tw 3 + u 4 )/s 7 , (17) 

|M^ X J 2 = 4096(145408m 10 - 1024mjjf - 1024m*u - 2368m^ 2 - 6400m^w - 2368m^ 2 
+ 16m 4 t 3 - 208m 4 t 2 w - 208m 4 tu 2 + 16m 4 w 3 + 24m 2 t 4 + 72m 2 c t 3 u + 96m 2 tV 
+72m 2 c tu 3 + 24m 2 w 4 - t 4 u - 3t 3 u 2 - 3t 2 u 3 - tu 4 ) / (3s 7 m 2 ). (18) 

For the known spin-triplet P-wave states Xcj(lP)(J = 0, 1, 2) we find 

cx(e + + e- - J/V(15) + Xco(lP)[Xci(lP), Xc 2 (lP)]) = 6.7[1.1, 1.6] fb, (19) 

a(e+ + e" -> ^(25) + Xco(lP)[Xci(lP), Xca(lP)]) = 4.4[0.74, 1.1] fb, (20) 
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For the exited spin-triplet 2P states Xcj(2P)(J = 0, 1,2), which are to be searched for, 
we find 

a(e + + e~ - J/^lS) + Xco(2P)[x c i(2P), x c2 (2P)]) = 9.1(1.6,2.2] fb, (21) 

a ( e + + e ^ ^2S) + Xc0 (2P)[ Xcl (2P), x c2 (2P)]) = 5.9[1.0, 1.4] fb. (22) 

In Eqs. (19-22), we see that the cross sections for the excited 2P states % c j(2P) are 
somewhat larger than that for the corresponding IP states Xcj(lP) in the nonrelativistic 
limit. This numerical result is due to the fact that we have chosen the first derivative 
of the wave function at the origin for the 2P states to be larger than that for the IP 
states, and actually the former could be slightly smaller than the latter, depending on the 
potentials that are used (see the values in the QCD (BT) model and other models in Ref.Q). 
Furthermore, another important effect comes from the relativistic corrections, which may 
lower the cross sections for the %cj(2P) states. E.g., if we take the charm quark mass to 
be 2m c = M(2P) « 4 GeV for the 2P states, then the cross sections will be substantially 
lower. Despite of these uncertainties, we expect that the cross sections for the 2P states 
should be comparable to that for the IP states. 



C. The r] c (nS)h c (mP) production 

To search for the spin singlet P-wave charmonium h c is certainly interesting. Recently, 
CLEO has found evidence for h c in the ip(2S) — > 7i°h c decay followed by h c — > ^rjc with a 

lyperfine splitting of about 1.0 MeV measured 



2 1[ |. It will also be interesting to search for the 



mass o£M(h c ) = 3524.4 ±0.6 ±0.4 MeV and 
in both the 7] c exclusive and inclusive analysis 
h c in e + e~ annihilation at y/s = 10.6 GeV in the recoil spectra of charge parity C=±l states 
such as T] c , XcO) Xcii Xc2- F° r the production of one spin-singlet S-wave state T] c (nS) and one 
spin-singlet P-wave state h c (mP), differing from [9j, we find |M| 2 to be not vanishing but 
given by Eq. (123]) : 

\M Vchc \ 2 = 2048(8192m^ ± 1024m^± 1024m^± 320m^ 2 ±128m^u± 320m^ 2 

- 16m^ 3 - 112m^ 2 u - 112m^u 2 - 16m^ 3 ± 2m 2 t 4 - 4m 2 c t 3 u - 12m 2 tV 

- 4m 2 /;w 3 + 2m 2 w 4 - t 4 u - 3tV - 3tV - tw 4 )/((8m 2 - t - u) 7 m 2 c ). (23) 
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For the spin-singlet P-wave states h c (lP) and h c (2P) we find 



a(e + + e~ -> ^(15) + /i c (1P)[/i c (2P)]) = 0.73[0.99] fb, (24) 

a(e + + e" -> r] c (25) + /i c (1P)[/i c (2P)]) = 0.48(0.65] fb, (25) 
a(e + + e" -> ?] c (35) + h c (lP)[/i c (2P)]) = 0.41(0.56] fb. (26) 

D. The Xcj(nP)(J = 0, 1, 2)h c (mP) production 

In the two P-wave (nP and mP) case, the cross section for e + + e~ — > 7* — > Hi(nP) + 
H 2 (mP) is 

<r(e+(pi) + e-(p 2 ) -> if^) + # 2 (p 4 )) = 
27ra 2 a 5 2 |P;(0)(V^-2 S K + ml) + (mI-m^ fi 

9^ L |M| dcos *- (27) 

For the Xcj^c production |M| 2 reads 

\M Xc0h f = 1638A(t 2 + u 2 -32m 4 c )(16m 2 c -3t-3u) 2 /(3(8m 2 c -t-u) 7 m 2 c ), (28) 



M Xclh f = -4096(8192m*° - 3072m** - 3072m s c u + 2048m^ 2 + 3584m^u 

+ 2048mj!u 2 - 16m 4 t 3 + 144m 4 t 2 w + 144m 4 tw 2 - 16m 4 w 3 - 52m 2 t 4 

- 128m 2 t 3 w - 152m 2 tV - 128m 2 tw 3 - 52m 2 w 4 + t A u 

+ 3tV + 3tV + tu 4 )/((8m 2 c -t-u) 7 m 4 c ), (29) 



\ M xc2hf = -8192(544m 4 + 72m 2 t + 72m 2 w + 3t 2 + 6tu + 3m 2 ) 

x (32m 4 c -t 2 -u 2 )/(3(8m 2 c -t-u) 7 m 2 c ), (30) 

and the corresponding cross sections are 

<r(e+ + e -> Xc o(lP) + K{1P) [h c (2P)\) = 0.22(0.31] fb, (31) 

a(e + + e" - Xcl (lP) + /i c (lP)[/i c (2P)]) = 1.0(1.4] fb, (32) 

(7(6+ + e~ -> x c2 (lP) + /i c (1P)[/i c (2P)]) = 0.063(0.085] fb. (33) 
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E. The tp(nS) 1 D 2 {mP) production 



The cross section for e + + e — > 7* — > Hi(nS) + H 2 (mD) process is formulated as 



(j(e + (pi) + e (p 2 ) -> ^(pg) + if 2 (p 4 )) = 

27r« 2 «2 1 i? 5 (0)| 2 |^(0) I V s2 - 2s(m§ + ml) + (m 2 - m 2 ) 2 



27m 2 s 2 



/: 



M| 2 dcos#, (34) 



I M J/ffl2 I = -327680(16m 2 - s) 2 (512m^ - 32m 4 s - 128m 4 t - 128m 4 w 

+ 8m 2 st + 8m 2 sw + 8m 2 t 2 + lQm 2 c tu + 8m 2 c u 2 - st 2 - su 2 ) / (3m 2 s 8 ) , (35) 

and the numerical result is 

a(e + + e~ -> 7/^(15) [2s] + 1 £> 2 ) = 0.19[0.12] fb. (36) 

F. The X cj{nP){J = 0, 1, 2) 3 D'j(mD)(J' = 1,2,3) production 

The cross section for e + + e~ — > 7* — > Hi(nP) + H 2 (mD) process is formulated as 



a(e + ( Pl ) + e-(p 2 ) - ifr fo,) + if 2 (p 4 )) = 

27ra 2 a 2 | J R P (0)| 2 |^(0)| 2 v /s 2 - 2s(m§ + m 2 ) + (m 2 - m 2 ) 2 '' J 



27 mis 



2 



M\ 2 dcos9, (37) 



1 



M Xcl 3 Dl I = 20480(113393664mJ 2 -49483776m^ s- 28348416m^- 28348416mJ M 

+ 7436288m^s 2 + 12601344m®st + 12601344m|5sw + 1829376m^ 2 

+ 3428352mj?fu + 1829376m^ 2 - 413120m^s 3 - 912960m^t 2 

- 1932032m°s 2 t - 1932032m^s 2 M - 1324416m^u - 912960m^sM 2 

+ 5600m 4 s 4 + 114816m 4 s 3 t + 114816m 4 s 3 w + 165024m 4 s 2 t 2 + 152960m 4 s 2 tw 

+ 165024m 4 sV + 144m 2 s 5 - 2312m 2 s 4 t - 2312m 2 s 4 w - 12184m 2 s 3 t 2 

- 4336m 2 s 3 tw - 12184m 2 sV + 289s 4 t 2 + 289sV)/(5m 4 s 9 ), (38) 
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I M Xcl iD2 | = 2560 (58589184m \ A - 2293760m^ 2 s - 14647296m 1 : 2 t - 4286464m 10 s 2 

- 14647296m 42 w + 6922240m^°st + 6922240m 40 sw + 2502656m^°t 2 

- 1343488m^°tM + 2502656m 4 V + 576512m®s 3 - 1119744m*s 2 i 

- 938496m®st 2 - 1119744m®s 2 w + 146432m^u - 938496m® su 2 

+ 4096m^ 4 + 72192m^ 3 t - 1248m 4 s 5 - 5984m 4 s 3 t 2 - 6080m 4 s 3 tw 
+ 72192m°s 3 w + 224m 2 s 4 t 2 + 119488m^s 2 t 2 + 40960m^s 2 tM 

- 1792m 4 s 4 t + 119488m°sV - 1792m 4 s 4 w - 5984m 4 sV 

- 48m 2 s 6 + 224m 2 sV + 3s 7 - 3s 5 t 2 + Qshu - 3s 5 u 2 )/(m 6 c s 9 ), (39) 

| M Xcl s Ds | = 40960(7446528mJ° - 756224m®s - 1861632m*t + 393856m^st 

- 1861632m®M - 62208m° s 2 + 393856m^w + 167552m^ 2 + 130304m^ M 
+ 167552m^ 2 + 9632m 4 s 3 - 27968m 4 s 2 t - 27968m 4 s 2 w - 39568m 4 st 2 

- 19328m 4 stw - 39568m 4 sw 2 - 288m 2 s 4 + 664m 2 s 3 t + 664m 2 s 3 w 
+ 752m 2 s 2 tw + 3120m 2 s 2 t 2 + 3120m 2 sV - 83s 3 1 2 

- 83sV)(16m 2 - s)/(5m 4 s 9 ), (40) 

M Xc2 3Dl | = 640(29807345664mJ 4 - 6599344128mJ 2 s - 7451836416mJ 2 t 

- 7451836416mJ 2 M + 182222848mJ°s 2 + 2217541632m* st + 2217541632m 10 
+ 647626752m 10 tw + 607666176m 10 t 2 + 607666176m*V + 43309056m^ 3 

- 233439232m®s 2 t - 233439232m®s 2 M - 205805568m^t 2 - 142774272mj!stu 

- 205805568m®sM 2 - 3057664m^s 4 + 10221568m|!s 3 t + 10221568m^ 3 w 
+ 24313856m^s 2 t 2 + 9732096m|!s 2 tM + 24313856m^V + 66368m 4 s 5 

- 163072m 4 s 4 t - 163072m 4 s 4 w - 1174848m 4 s 3 t 2 - 400m 2 s 6 

- 205696m 4 s 3 te - 1174848m 4 sV + 20768m 2 sV 

- 768m 2 s 4 tw + 20768m 2 sV + s 7 - sH 2 + 2sHu - s 5 u 2 )/(15m 6 c s 9 ), (41) 
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M Xc2 a D2 I = 81920(13979648m 42 - 1540096m 40 s - 3494912m 40 t - 3494912m 40 w 

- 62400m^s 2 + 724224m*st + 724224m*su + 303232m^ 2 + 267264mj?*u 
+ 303232m^ 2 + 14208m^s 3 - 55440m^s 2 t - 55U0m b c s 2 u - 71376m^t 2 

- 38304m°stu - 71376m^sM 2 - 624m 4 s 4 + 1872m 4 s 3 t + 1872m 4 s 3 w 
+ 6018m 4 s 2 t 2 + 1824m 4 s 2 tw + 6018m 4 sV + 9m 2 s 5 - 24m 2 s 4 t 

- 24m 2 s 4 M - 219m 2 s 3 t 2 - 30m 2 c sHu - 219m 2 s 3 w 2 

+ 3s 4 t 2 + 3sV)/(m 4 s 9 ), (42) 

| M Xc2 3 Ds | = 5120(883687424m 44 - 77856768m 42 s - 220921856m 42 t 

- 220921856m 42 w - 2650112m 40 s 2 + 40353792m 40 st + 40353792m 40 sw 
+ 19030016m 40 t 2 + 17170432m 40 tw + 19030016m*V + 1043456m^s 3 



- 2798592m*s 2 t - 2798592m^s 2 w - 3923968m^st 2 - 2240512m*stu 

- 3923968m^sM 2 - 95744m°s 4 + 89088m° s 3 t + 89088m°s 3 w 
+ 332416m^s 2 t 2 + 3A816m 6 c sHu + 332416m°sV + 3968m 4 s 5 

- 1152m 4 s 4 t - 1152mVw - 12928m 4 s 3 t 2 + 3584m 4 s 3 tw 



12928m 4 sV - 80m 2 s 6 + 208m 2 s 4 t 2 - 128m 2 s 4 tw + 208m 2 s 4 w 2 



+ s 7 - sH 1 + 2sHu - sV)/(5m°s 9 ). (43) 

The numerical results are listed as follows, where 5\ = 3 Di, 5 2 = 3 D 2) 5% = 3 _D 3 , 

<y{e + + e -> X ci\Xc2\ + S ± ) = 0.080(0.041] fb, (44) 

<y{e + + e -> Xci[Xc 2 ] + 5 2 ) = 0.099[0.084] fb, (45) 

<y{e + + e" -> Xci[x c2 ] + 5 3 ) = 0.041[0.0099] fb, (46) 

As a summary of the above results, we show the differential cross sections (the angular 
distribution functions) for different double charmonium production processes at leading 
order in NRQCD (QED contributions are not included) in Table I, and the corresponding 
graphs in Fig. 1-7. 
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III. QED PROCESSES INCLUDING J PC = 1 STATES 

Since the J PC = 1 cc can be produced via a single photon, the QED contribution 
may be significant or even comparable to the QCD contribution in some exclusive processes 
involving one J PC& = 1 charmonium state. These QED effects are considered in j^], H ]. 
There are six Feynman diagrams for the QED process at a 4 order, and only two represent 
7* — > ccy* — > (cc)i — cc, which is dominant and has been calculated in 0|. In this paper 
we include all the eight diagrams to get the full result at order a 4 , though the contributions 
of other six diagrams are numerically small. Using the notation in section II, we re-express 
below the analytical formulas of the exclusive processes including both QCD and QED 
contributions. 

A. The ip(nS)r] c (mS) production 

The cross section for e + + e~ — > 7* — > Hi(nS) + H 2 (mS) is now changed to 
a(e + ( Pl ) + e-{p 2 ) -> H^ps) + H 2 (p A )) = 

a*\R ns (0)\z\R ms (0)\Ws 2 - ggg + gD + K ~ [' ,^,2 , (A7 , 

288m2vr S 2 JJM\ ax, ^t) 

where x = cos 9 and 9 is the scattering angle between p± and p%. And |M| 2 is 

- 2 2048(s- 16m 2 c ){16am 2 c + A8a s mc 2 + 3as) 2 
lMl ~ 81^ (48) 

The numerical results become 

a(e + + e~ -> J/^(1S) + r/ c (lS) [r] c (2S), ^(35)]) = 6.6(4.3,3.7] fb, (49) 

a(e+ + e" -> tf;(2S) + rj c (lS) [ Vc (2S),r] c (3S)}) = 4.3[2.8, 2.4] fb, (50) 

where the values in the brackets are the cross sections respectively for rj c (2S) and rj c (3S) 
with recoiled ip mesons. 
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B. The ifi(ns) + Xcj(mp) Production 



The cross section for e + + e — > 7* — > H\{nS) + H2(mS) is changed to 
a{e + { Pl ) + e-(pa) -> #i(p 3 ) + H 2 (p 4 )) = 

j \M\ 2 dx. (51) 



a 2 \R ns (0)\ 2 \R' mp (0)\Ws 2 - 2s(m 2 + m|) + (m 2 - ^ f\ 



96m 2 irs 2 

And now |M| 2 for ipXco,4'Xci,4'Xc2 production, are given in Eq. ( l52l . Eq. ( 153]) . Eq. ([5 
respectively. 

(M^J 2 = (2048tt 2 ((589824(x 2 - l)mf - 4096s(73x 2 + 25)m 8 c - 2048s 2 (x 2 - l)m 6 c 
+128s 3 (25x 2 + 17)m 4 c + 16s 4 (x 2 - l)m 2 - 9s 5 (x 2 + l))a 2 
+48a s m 2 (73728(x 2 - l)m 8 c - 3584s(13x 2 + l)mj! + 64s 2 (71x 2 + 55)m 4 
+56s 3 (5x 2 + l)m 2 - s 4 (25a; 2 + 29))a + 144a 2 m 2 (36864(x 2 - l)m 8 c 
-256s(109x 2 - ll)77ijj + 128s 2 (41x 2 + 22)m 4 - 4s 3 (61x 2 + 101)m 2 
+s\x 2 - l))))/(243m 4 s 6 ). (52) 



\M^ Xcl | 2 = (40967r(192aa,(9216(x 2 - l)m° - 64s(31x 2 + 19)m 4 + 8s 2 (9x 2 - 19)m 2 
+s 3 (x 2 + 25))m 4 + 2304a 2 (1152(x 2 - l)m 6 c - 8s(49x 2 + l)m 4 
+4s 2 (9a; 2 + 5)m 2 - s 3 (a; 2 + l))m 4 + a 2 (294912(a; 2 - \)mf 
-2048s(13s 2 + 37)m* - 19456s 2 m^ + 32s 3 (x 2 - 2Z)m% + 144sVm 2 
-9s 5 (x 2 + l))))/(81m 4 s 6 ). (53) 



\Mf x J 2 = (4096tt 2 ((1474560(x 2 - l)m* u + 2048s(31x 2 - 233)m° - 1024s 2 (32a; 2 + l)m b c 
-32s 3 (35x 2 - 101)m 4 + 16s 4 (19x 2 - 10) - 9s 5 (x 2 + l))a 2 
+48a s m 2 ( 184320 (x 2 - l)m 8 c - 256s(59x 2 + U3)m 6 c - 32s 2 (59a; 2 - 89)m 4 
+4s 3 (61a; 2 + 29)m 2 - s 4 (7x 2 + ll))a + 144a 2 m 2 (92160(:r 2 - l)m 8 c 
-128s(149x 2 + 53)m^ + 64s 2 (25x 2 + 23)m 4 - 32s 3 (2x 2 + l)m 2 
+s 4 (x 2 - l))))/(243m 4 s 6 ). (54) 

The numerical results of spin-triplet P-wave states X cjO-P)(J — 0, 1, 2) are 

a(e + + e- -> J/^IS) + Xc0 (lP)[Xci(lP),Xc2(lP)}) = 6.9[1.0, 1.8] fb, (55) 
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a(e+ + e - if,(2S) + Xc0 (lP)[ Xc i(lP), X c2 (lP)]) = 4.5[0.7, 1.1] fb, (56) 
And for the exited spin-triplet 2P states Xcj(2P)(J = 0, 1, 2), the results turn to be 

a(e + + e~ - 7/^(15) + Xc o(2P)[x c i(2P), * c2 (2P)]) = 9.4[1.4, 2.4] fb, (57) 

a(e + + e" - ^(25) + x c0 (2P)[ Xcl (2P), Xc2 (2P)]) = 6.2[0.9, 1.6] fb. (58) 

C. The ip(ns) + 1D2 Production 

The cross section for e + + e~ — > 7* — > H-y{nS) + H 2 (mD) process is 
^(e + (pi) + e~(p 2 ) -> ffi(p 3 ) + P2G04)) = 

5 Q 2 |P S (0)| 2 |P£(0)| Vg - 2s(mj + m 2 ) + (mg ~ ^!) 2 1^,2 . ,, Q , 

192/^2 J VVl\ax, [M) 



and 



|2 40967T 2 (s-16m2) 3 (32am2 + 96a s m2 + 3sa) 2 (x 2 + l) 
Mj/^id 2 I = , (60) 



c 



and the numerical result becomes 

a(e + + e" -> 7/^(15) [2s] + 1 P> 2 ) = 0.21[0.13] fb. (61) 



D. The rj c (ns) + ip3Di Production 

The cross section for e + + e~ — > 7* — > Hi(nS) + H 2 (mD) process is formulated as 



cr(e + (pi) + e-(p 2 ) -> ifi(p 3 ) + P 2 (p 4 )) = 
5a 2 |P 5 (0)| 2 |P£(0)| V* 2 - 2s(m§ + m 2 ) + (m 2 - m 2 ) 2 



192m27rs 2 



c 



y |M| 2 dx, (62) 



and 



|M„ c 3 Dl | 2 = 512tt 2 (s - 16m 2 )(48a s (64m 2 - 9s)m 2 + a(102Am 4 c - UAsm 2 c - 15s 2 )) 2 

(x 2 + l)/(1215m^s 6 ), (63) 

and the numerical result becomes 

a(e + + e- ^ rj c (lS) + 3 D ± ) = 0.17 fb. (64) 
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E. conclusion 

In this section we have considered the QED contribution to the exclusive processes in- 
volving one J pc = 1 charmonium state (e.g., ip(nS) and ijizDi)- We find that in general 
by adding QED contribution the cross section can be changed by the order of ten percent. 
We list the results below with y/s = 10.6GeV, a s = 0.26, m c = l.bGeV. The cross sections 
of i[}(ns) + rj c (ms) are increased by 20 percent, that of ip(ns) + Xcj are increased by 4, -5, 8 
percent for J = 0,1,2 respectively, that of %fj{ns) + 1 D 2 is increased by 11 percent, and 
that of rj c (ns) + 3 D\ is increased by 9 percent. In |9|] the authors also considered the QED 
process but only the two dominant diagrams were included. Our results are in agreement 
with their results in most processes when we choose the same parameters as theirs, except 
for the process of r] c + ipi^Di) production, for which they obtained 41% enhancement with 
QED effects whereas we get 19% with the same parameters. Here our analytical expression 
also differs from theirs. 

We show the differential cross sections (the angular distribution functions) for different 
double charmonium production processes including both QCD and QED contributions in 
Table II, and the corresponding graphs in Fig. 8. Ratios of production cross sections of 
various double charmonia to that of J/ip + r\ c in e + e~ annihilation at y/s = 10.6 GeV are 
listed in Table III. 



IV. DISCUSSION AND SUMMARY 

In this paper, we make predictions for various double charmonia production processes in 
e + e _ annihilation at y/s = 10.6 GeV with B factories, based on a complete leading order 
calculation including both QCD and QED contributions. In particular, we aim at searching 
for excited charmonium states in these processes. From the obtained results we make the 
following observations: 

The calculated relative production rates for e + e~ — > jib{2S)r] c , ip(2S)x c o,i J (^S)r] c (2S) are 



roughly compatible with the new Belle measurements [8| (see also j^]), assuming the decay 
branching ratios into charged tracks are comparable for r] c ,T] c (2S), and Xco- 

The calculated relative production rates for e + e~ — > J j ' ipXcoi^P) , J/ipi]c(3S) and e + e _ — > 
tp(2S)Xco{2P),'4 ) (2S)Vc(3S) ar e large, and these two states may be observable in the mass 
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range m 2 p = 3.90 — 4.00 GeV and m 3 s = 3.95 — 4.10 GeV respectively (see, e.g. [6j). 
Both of them are above the OZI (Okubo-Zweig-Iizuka) allowed thresholds, but unlike the 
Xco(2-P), the r/ c (3S) can not decay to DD pair, which may distinguish between these two 
states experimentally. 

The calculated relative production rate for e + e~ — > h c r) c is not zero. This differs from the 
result given in 9j, but agrees with |22j. Hence the h c might be observable via this channel 
with high statistics in the future. Aside from this process, e + e~ h c \ c \ is also very hopeful 
in finding the h c meson, since the Xci has a large branching ratio (larger than 30%) decaying 
into ■yJ/tp, and the J /ip can be easily detected by the n + n~ signal. Moreover, the calculated 
cross section for e + e~ — > h c x c i is about 1.0 fb, not very small and much larger than that for 
h c Xco an d h c Xc2 production rates. 

We show the differential cross sections (the angular distribution functions) for different 
double charmonium production processes in Table I (not including QED contribution) and 
Table II (including both QCD and QED contributions), and the corresponding graphs in 
Fig. 1-7, and Fig. 8. 



As a whole, our results a gree with those in 9] and 22J (also agree with our previous 



result for ■//<.'//..( \ ,.n. Xci ? Xc2)[lfll)) if using the same parameters. For the QED part, in 
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22| only two dominant diagrams are taken into account, while in this paper all diagrams 



are considered, but the numerical contributions of the remaining four diagrams are small. 
However, there still exists a difference in the result for the e + e~ — ► r] c + ip( s Di) production. 
As for numerical results, since we use a larger value of the strong couplingconstant a s = 0.26 
(corresponding to fi = 2m c = 3.0 GeV and A^-=338 MeV) than that in 9, 22], our predicted 
cross sections are in general larger than that in j^], [22]. Moreover, we use the charmonium 
wave functions and their derivatives at the origin (including the ground state and excited 
states) from the BT potential model calculation [20j but not from the experimental values 
of leptonic decay widths, etc., as in 0, Q|, and this may further enlarge our predicted cross 
sections. These parameters may not be the best choice, but, at present, since we do not 
have enough available data for higher excited charmonium states, using the potential model 
calculation may still be a reasonable and tentative choice. We view these as theoretical 
uncertainties in our approach for the leading order calculations. 

As emphasized above, this paper aims at searching for excited charmonium states in e + e~ 
annihilation at t/s = 10.6 GeV with B factories. After the main part of the results were 
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presented in [23J (with a more complete leading order calculation including both QCD and 
QED contributions for various double charmonium production processes being added in its 
present form), a number of new experimental and theoretical results have appeared recently. 

1. The double charmonium production in e + e~ annihilation at B factories has been 
confirmed by the BaBar Collaboration (see Ref. 24|) with comparable cross sections to that 
observed by Belle. Theoretically, the large gap between experiment and theory could be 



largely narrowed by the next to leading order QCD radiative corrections 
corrections (see, Refs. 



26 



251 ] and relativistic 



27] and references therein) in the framework of nonrelativistic 

n 

QCD, and other possible approaches (see, e.g. [28]). 

2. Belle has observed the X(3940), a new charmonium state or charmonium-like state, 
in e + e- -> J/^X(3940) with M x = (3.940 ± 0.012) GeV|2a|. While its main decay mode 
is X(3940) — > DD*, no signal is found for X(3940) — > DD. This rules out the possibility of 
X(3940) being the x c o(2P) state. Furthermore, the X(3940) is unlikely to be the Xd(2P) 
state, since no signal for the XdO-P) is found. This is in line with our calculation, which 
shows that the cross section for e + e~ — > J/i^XciO-P) is much smaller than e + e _ — > J/i/jt) c . 
Finally, the X(3940) could be the r] c (3S) state. According to our calculation (see Eq.(13)), 
the production cross sections for i] c (lS), r) c (2S),r) c (3S) are respectively 5.5, 3.6 ,3.1 fb, and 
the relative rates are roughly consistent with the signal yields N = 235±26, 164±30, 149±33 
events for r) c (lS), rj c (2S), X(3940) {29]]. This might be viewed as a support to interpreting 
the X(3940) as the r] c (3S). However, the remaining problem is how to understand its low 
mass of X(3940) if it is the r) c (3S), which is lower than potential model predictions by 50- 
120 MeV. But this could be explained by the coupled channel effects that the coupling of 
r] c (3S) to the + and 0~ charmed meson pair (in S-wave) will lower the mass of i] c (3S) 30]. 

3. Belle has not observed the X(3872) in e + e~ — ► J/^X(3872). This implies that the 
X(3872) is unlikely to be a conventional 0~ + or ++ charmonium, since in our calculation 
they may have relatively large rates to be observed in double charmonium production. On 
the other hand, a 1 ++ or 2 h charmonium could be possible for X(3872), since they have 
relatively small production rates. Of course, the nature of X(3872) needs clarifying by other 
more relevant experiments (see, e.g. 5j), aside from the e + e~ annihilation processes. 

4. Belle has very recently observed a new state, the X(4160), in the process of double 



charm production e + e" -> J /if) + X(4160) followed by X(4160) -> £>*£>* |3lJ. Possible 



interpretations for the X(4160) are discussed in 32j with emphasized possible assignments 
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of 2 1 D 2 , Tj c (AS), Xco(3-P) charmonium states and related problems. 

In conclusion, we find that the double charmonium production processes in e + e~ anni- 
hilation at B factories are very useful tools in searching for excited charmonium states or 
charmonium- like states. The leading order NRQCD calculation might hopefully provide a 
useful guide for the relative production rates, but not the absolute rates themselves. A 
systematical study for the QCD radiative corrections and relativistic corrections for dif- 
ferent processes are apparently needed. On the other hand, studies of charmonium spec- 
troscopy including charmonium masses, decays, and coupled channel effects, and the new 
type charmonium-like states are also very desirable. 

Acknowledgments 

We thank P. Pakhlov for discussions on the Belle data and possible implications, and E. 
Braaten and J. Lee for communications on some calculations. This work was supported in 
part by the National Natural Science Foundation of China (No. 10421503, No. 10675003), 
the Key Grant Project of Chinese Ministry of Education (No. 305001), and the Research 
Found for Doctorial Program of Higher Education of China. 



T. Skwarnicki, Int.J.Mod.Phys. A19, 1030(2004), hep-ph/0311243 



S.L. Olsen, Int.J.Mod.Phys. A20, 240(2005), hep-ex/0407033, 



C. Quigg, hep-ph/0403187; Nucl.Phys.Proc.Suppl. 142, 87(2005). 



F.E. Close and P.R. Page, Phys.Lett. B578, 119(2004), hep-ph/0309253 



E. Swanson, hep-ph/060lTTol 

S. Godfrey, N. Isgur, Phys. Rev. D32, 189 (1985). 

Belle Collaboration, K. Abe et al, Phys. Rev. Lett. 89, 142001 (2002). 



Belle Collaboration, K. Abe et al, hep-ex/0407009, Phys. Rev. D70, 071102 (2004) 



E. Braaten and J. Lee, Phys. Rev. D67, 054007 (2003). 

K.Y. Liu, Z.G. He and K.T. Chao, Phys. Lett. B557, 45 (2003). 

K. Hagiwara, E. Kou, C.-F. Qiao, Phys. Lett. B570, 39 (2003). 



18 



[12] G. T. Bodwin, J. Lee, E. Braaten, Phys. Rev. D67, 054023 (2003); Phys. Rev. Lett. 90, 
162001 (2003); 

A. B. Kaidalov, JETP Lett. 77, 349 (2003), Pisma Zh. Eksp. Teor. Fiz. 77, 417 (2003), 



hep-ph/0301246 



A. V. Berezhnoy, A. K. Likhoded, Phys. Atom. Nucl. 67, 757(2004) ; Yad.Fiz. 67, 778(2004), 
|hep-ph/0303145| 

K.-Y. Liu, Z.-G. He, K.-T. Chao, Phys. Rev. D68, 031501(R) (2003); 

A. V. Luchinsky, Phys. Atom. Nucl. 68, 1400(2005); Yad. Fiz. 68, 1456(2005), 



hep-ph/0305253; 



S. J. Brodsky, A. S. Goldhaber, J. Lee, Phys. Rev. Lett. 91, 112001 (2003); 
B. L. Ioffe, D. E. Kharzeev, Phys. Rev. D69, 014016(2004), |hep-ph/0306062~t 



S. Fleming, A. K. Leibovich, T. Mehen, Phys. Rev. D68, 094011(2003) , |hep-ph/0306139; 

C.-F. Qiao, J.-X. Wang, Phys. Rev. D69, 014015(2004), |hep-ph/0308244| 

S. Dulat, K. Hagiwara and Z.H. Lin, Phys. Lett. B594, 118(2004), [hep-ph/04022"30 

K. Cheung and W.Y. Keung, Phys. Rev. D69, 094026 (2004). 



J.P Ma and Z.G. Si, Phys. Rev. D70, 074007 (2004), piep^ph/0405111 



[13] P. Cho and A.K. Leibovich, Phys. Rev. D 54, 6690 (1996). 

[14] F. Yuan, C.F. Qiao, and K.T. Chao, Phys. Rev. D 56, 321 (1997); Phys. Rev. D 56, 1663 
(1997). 

[15] S. Baek, P. Ko, J. Lee, and H.S. Song, J. Korean Phys. Soc. 33, 97 (1998). 

[16] KY. Liu, Z.G. He, and K.T. Chao, Phys. Rev. D69, 094027 (2004). 

[17] G.T. Bodwin, E. Braaten, and G.P. Lepage, Phys. Rev. D51, 1125 (1995). 

[18] P. Cho and A.K. Leibovich, Phys. Rev. D53, 150 (1996); 53, 6203 (1996); P. Ko, J. Lee and 

H.S. Song, Phys. Rev. D54, 4312(1996); ibid D60, 119902(1999). 
[19] J.H. Kiihn, J. Kaolan and E.G.O. Safiani, Nucl. Phys. B157, 125 (1979); B. Guberina, J.H. 

Ktihn, R.D. Peccei and R. Riickl, Nucl. Phys. B174, 317(1980). 
[20] E.J. Eichten and C. Quigg, Phys. Rev. D52, 1726 (1995). 

[21] P. Rubin et al.(CLEO Collaboration), Phys.Rev. D72, 092004(2005), |hep-ex/0508037 

[22] E. Braaten and J. Lee, Phys. Rev. D72, 099901(E) (2005). 

[23] K.Y.Liu, Z.G.He and K.T.Chao, | arxiv:hep-ph/0408141| y2. 

[24] B. Aubert et al, The BABAR Collaboration, Phys.Rev. D72, 031101(2005). 



19 



[25] Y.J. Zhang, Y.J. Gao and K.T. Chao, Phys. Rev. Lett. 96, 092001(2006); see also Y.J. Zhang 
and K.T. Chao, Phys. Rev. Lett.98, 092003 (2007). 



[26] G.T. Bodwin, D. Kang and J. Lee, Phys. Rev.D74, 014014 (2006)(hep-ph/0603186 ); 



hep-ph/0603185; G.T. Bodwin, D. Kang, D. Kim, J. Lee, and C. Yu, hep-ph/0611002 



[27] Z.G. He, Y. Fan, and K.T. Chao, Phys. Rev. D75, 074011 (2007). 

[28] J.P Ma and Z.G. Si, |arXiv:hep^ph70 608221 ; V.V. Braguta, A.K. Likhoded, and A.V. Luchin- 
sky, Phys. Rev. D72, 074019 (2005); [ailQv:hep-ph/0602047l |arXiv:hep-ph/06lT02lj A.E. Bon- 
dar and V.L. Chernyak, Phys. Lett. B612, 215 (2005); D. Ebert and A.P Martynenko, Phys. 
Rev. D74, 054008 (2006); H.-M. Choi and C.-R. Ji, arXiv: 0707.1173 [hep-ph]. 



[29] P. Pakhlov (Belle Collaboration), hep-ex/0412041 



[30] E. Eichten, K. Lane, and C. Quigg, Phys. Rev. D73, 014014 (2006); 73, 079903(E) (2006). 
[31] I. Adachi et al. (Belle Collaboration). !arXiv:0708.3812l [hep-ex]. 
[32] K.T. Chao, arXiv:0707.3982l hep-ph] . 



20 



TABLE I: Differential cross sections-angular distribution functions for double charmonium pro- 
duction in e + e~ annihilation at ^/s = 10.6 GeV (QED contributions not included; see text for the 
input parameters). 



Differential cross section 


Angular distribution function (fb) 


da(e + + e' -> J lib + riAlS)) I dcosO 


2.06(1 + cos 2 #) 


da(e + + e~ -► JM + ri c (2S))/dcos9 


1.34(1 + cos 2 6») 


da(e + + e~ -► JM + r] c (3S))/dcos9 


1.16(1 + cos 2 6») 


dcr(e + + e~ — > JM + Y r n)/dcos9 


3.09(1 + 0.252cos 2 fl) 


dcr(e + + e~ ^ J lib + \ c i)/dcos9 


0.457(1 +O.698cos 2 0) 


da(e + + e~ — ► J 'M + Y C 2)/dcos0 

\ IT' /V "--^ / / 


0.870(1 - 0.198cos 2 6l) 


afo"(e + + e~ — > ri r + hc(lp)) / dcos9 


0.480(1 - 0.726cos 2 6>) 


da(e + + e~ — ► ri c + hc(2p)) / 'dcosO 


0.653(1 - 0.726cos 2 6») 


da(e + + e~ — ► y c i + hc(lp)) / dcos9 

\ /V.L--L \ 1 / / 1 


0.483(1 +O.19Ocos 2 0) 


da(e + + e~ — ► y c i + hc(2p)) 1 dcosO 

\ /V.L.J- \ Jr / / / 


0.657(1 +0.190cos 2 #) 


da(e + + e~ — > y C 2 + hc(lp)) I dcosO 


0.0235(1 + cos 2 6>) 


da(e + + e~ — > y C 2 + hc(2p)) 1 dcosO 


0.0306(1 + cos 2 6>) 


a(e + + e - ->■ J/^(1S) + 1 D 2 )/dcos9 


0.0694(1 + cos 2 #) 


da{e + + e~ — > % c i + 5\)/dcos9 


0.0307(1 + 0.922cos 2 #) 


da{e + + e~ — > Xci + 8z)/dxx)s6 


0.0545(1 - 0.273cos 2 #) 


da{e + + e~ — > Xci + 5s)/dcos9 


0.0163(1 + 0.760cos 2 6») 


da{e + + e~ — > Xc2 + 5\)/dcos9 


0.0161(1 +0.952cos 2 #) 


da(e + + e" -> Xc2 + 5 2 )/dcos9 


0.00331(1 + 0.799cos 2 6>) 


da(e + + e~ — > Xc2 + S 3 )/dcos9 


0.00530(1 - 0.188cos 2 #) 
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TABLE II: Differential cross sections-angular distribution functions for double charmonium pro- 
duction in e + e~ annihilation at ^/s = 10.6 GeV with both QCD and QED contributions (see text 
for the input parameters). 



Differential cross section 


Angular distribution function (fb) 


da(e + + e~ -> J lib + nJlS)) 1 dcosO 


2.47(1 + cos 2 #) 


da(e + + e" -► J/Y> + <q c {2S)) / dcosO 


1.62(1 + cos 2 6>) 


da(e + + e~ J/Y> + r? c (3S"))/(fcosf9 


1.39(1 + cos 2 <9) 


cfo"(e + + e~ — ► J/V' + Xco)/dcos9 


3.18(1 + 0.265cos 2 #) 


da(e + + e~ — > J/V> + Xci)/dcos9 


0.426(1 + 0.751cos 2 #) 


da(e + + e~ — > J/V' + Xci)/dcos9 


0.929(1 -0.161cos 2 f9) 


da(e + + e~ ->■ J/*b{lS) + 1 D 2 )/dcos6 


0.0770(1 + cos 2 6) 
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TABLE III: Ratios of production cross sections of various double charmonia to that of J ' /ip + r\ c in 
e + e~ annihilation at \fs = 10.6 GeV. 





n c (lS, 25, 35) 


Y c0 flP,2P) 


y c i(1P,2P) 


y c2 (1P,2P) 


h c (lP,2P) 


3 D! 


3 D 2 


3 ^3 




ib(lS) 


1.0,0.65,0.56 


1.05,1.4 


0.15,0.21 


0.27,0.36 










0.03 


ib(2S) 


0.65,0.42,0.36 


0.68,0.94 


0.11,0.14 


0.17,0.24 










0.02 


Vc(lS) 










0.11,0.15 


0.025 








r?c(2S) 










0.07,0.10 


0.016 








r? c (3S) 










0.06,0.08 










Xo(lP) 










0.03,0.05 










Xi(lP) 










0.15,0.21 


0.012 


0.015 


0.006 




^(IP) 










0.010,0.013 


0.006 


0.001 


0.0015 
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FIG. 1: Differential cross sections for e + e — > J/ifi + Xcj{f]c) 
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FIG. 2: Differential cross sections for e + e — > r/ c + h c 
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FIG. 3: Differential cross sections for e + e — > Xd + h c 
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FIG. 4: Differential cross sections for e + e — > Xc2 + h c 
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FIG. 5: Differential cross sections for e + e — > J/ip + 1 Z>2 
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FIG. 6: Differential cross sections for e + e — ► Xd + $J 
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FIG. 7: Differential cross sections for e + e — > \c2 + $j 
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